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SECTION  1 


INTRODUCTION 

This  report  details  the  coupling  of  the  CM3DT  Invlscld  nosetlp 
flow  field  code  to  the  STEIN  Invlscld  afterbody  flow  field  code.  The  CM3DT 
code,  described  In  References  1  and  2,  is  a  time-dependent  code,  in  which 
the  steady-state  solution  for  supersonic/hypersonic  flow  over  a  nosetlp  Is 
obtained  as  the  asymptotic  limit  of  the  unsteady  flow  problem.  The  STEIN 
code,  described  in  Reference  3,  Is  a  forward  marching  solution  of  the  steady 
supersonic  inviscid  flow  equations,  and  requires  the  specification  of  an 
initial  flow  field  data  plane,  which  must  be  obtained  from  a  nosetip  flow 
field  code. 

The  STEIN  code  was  specifically  developed  to  treat  complex  geometries, 
such  as  the  Space  Shuttle  Orbiter,  and  includes  the  ability  to  rigorously  treat 
embedded  shocks  that  occur  in  the  flow  field  surrounding  such  configurations 
(e.g.,  wing  shocks).  To  compute  these  complex  configurations,  STEIN  uses  a 
sequence  of  conformal  transformations  to  produce  a  body-fitted  computational 
coordinate  system. 

Geometries  in  the  STEIN  code  are  defined  using  the  QUICK  geometry 
system  (Reference  4),  which  allows  the  definition  of  complex  geometries  (such 
as  the  Shuttle  Orbiter)  through  specification  of  body  cross-sections  and 
meridional  profiles. 

The  goaT  of  this  effort  is  to  couple  the  CM3DT  nosetip  code  to  the 
STEIN  afterbody  code.  This  coupling  Involves  relating  the  CM3DT  and  STEIN 
coordinate  systems  and  developing  an  interpolation  procedure  for  defining 
the  Initial  data  required  for  the  STEIN  code  from  CM3DT  solutions.  In  addi¬ 
tion,  to  ensure  compatibility  between  CM3DT  and  STEIN,  the  QUICK  geometry 
system  has  been  Implemented  into  CM3DT,  allowing  a  single  definition  of  a 
vehicle  geometry  to  be  used  In  both  the  nosetip  and  afterbody  flow  field 
codes. 
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The  specific  goal  of  coupling  the  CM3DT  and  STEIN  codes  was  to 
obtain  Inviscld  calculations  on  the  Shuttle  Orblter  at  higher  angles  of 
attack  than  had  been  previously  possible.  Incorporation  of  the  QUICK 
geometry  package  In  CM3DT  simplifies  definition  of  the  Shuttle  Orblter 
nosetip  geometry,  using  existing  descriptions  of  the  Orblter  geometry  in 
the  QUICK  format. 

Section  2  of  this  report  describes  the  implementation  of  the 
QUICK  geometry  system  into  the  CM3DT  code  and  details  the  relationship  be¬ 
tween  the  respective  coordinate  systems.  The  generation  of  Initial  data 
for  the  STEIN  code  is  described  In  Section  3.  A  brief  guide  to  the  use 
of  the  QUICK/CM3DT/STEIN  system  of  codes  Is  provided  in  Section  4. 
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SECTION  2 


IMPLEMENTATION  OF  QUICK  GEOMETRY  PACKAGE 
INTO  CM30T  FLOW  FIELD  CODE 

The  QUICK  geometry  code  was  developed  at  Grumman  Aerospace  Corpora¬ 
tion  for  NASA  to  provide  a  means  of  defining  complex  vehicle  geometries  in  a 
relatively  simple  manner.  This  code  was  structured  to  allow  its  Incorporation 
into  flow  field  codes  without  modification.  Details  on  the  geometry  models 
used  in  QUICK  and  instructions  for  its  use  may  be  found  in  References  4  and  5. 
(Codes  to  simplify  the  required  Inputs  to  QUICK  ere  described  in  References 
6  and  7). 

The  QUICK  geometry  system  assumes  a  pitch  plane  of  geometric  sym¬ 
metry;  l.e.,  in  the  QUICK  Cartesian  coordinate  system  shown  In  Figure  2.1, 
the  Xg  -  Zg  plane  is  a  plane  of  symmetry.  The  geometric  Inputs  to  the  QUICK 
code,  which  will  not  be  discussed  in  this  report,  are  defined  In  the  (xq, 
yq ,  Zg)  coordinate  system. 

The  output  from  the  QUICK  code  used  for  geometric  definitions  in 
flow  field  codes  is  defined  in  a  cylindrical  coordinate  system.  Since  QUICK 
does  not  permit  the  body  radius  in  this  cylindrical  system  to  be  multi¬ 
valued,  provision  has  been  made  to  allow  for  a  "moving  axis",  which  serves 
as  the  center  of  the  (r,  0)  polar  frame  within  each  cross-section  Xg  = 
constant.  The  moving  axis  is  constrained  to  lie  in  the  xq  -  Zg  plane,  and 
is  defined  by  the  curve  zq  a  d  (xq ) ,  yq  -  0.  Figure  2.1  illustrates  the 
resulting  (r,  0,  Xg)  coordinate  system,  which  is  related  to  the  (xg,  yg,  Zg) 
Cartesian  system  through 

r  =  /yq2  +  (zQ  -  d)2  (2.1) 

0  *  tan'1  ZQ~  _  .  (2.2) 

*Q 

With  the  plane  of  symmetry  assumption  In  QUICK,  the  polar  angle  0 
is  constrained  to  the  range  -  */2  <_  0 
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In  the  CM3DT  nosetlp  flow  field  code  the  body  geometry  Is  defined  in 
a  cylindrical  coordinate  system  (x,  y,  <J>)»  which  is  related  to  a  corresponding 
Cartesian  system  (xc,  yc,  2C)  through 

xc  -  x  (2.3) 

yc  *  y  cos  $  (2.4) 

zc  -  y  sin  <i>  (2.5) 

as  shown  in  Figure  2.2;  the  polar  angle  $  is  constrained  to  the  range  0  <  $  £  2it. 

The  QUICK  and  CM3DT  coordinate  systems  are  related  by  the  specification 
of  the  parameters  #b0  and  *0,  which  are  input  to  the  CM3DT  code  as  XBQ  and  XHQ, 
respectively  (see  Section  4).  $b0  represents  the  Xq  value  of  the  most  forward 
nosetlp  geometry  point,  and  Is  a  second  value  of  xq  selected  by  the  user  of 
the  code  to  define  the  nosetip  computational  axis. 

The  CM3DT  computational  axis  is  defined  as  the  straight  line  that 
passes  through  the  points  xg  =  *b0«  *  0  Zq  “  d(*b0)  a  d-j  and  Xg  *  8Q,  y^  - 

0,  Zq  *  d  (XQ)  =  d2  in  the  QUICK  Cartesian  coordinate  system.  This  CM3DT  axis 
lies  entirely  within  the  xq  -  zq  plane  and  Is  oriented  at  an  angle  6^  relative 
to  the  xq  axis,  where 

.  drd, 

6  -  tan’1  ^-—2-  (2.6) 

1  xo  ’  xb0 

as  shown  in  Figure  2.3. 

This  angle  <5-|  is  similar  in  definition  to  the  angle  6^  (DELI)  used 
in  the  original  formulation  of  CM3DT,  as  described  in  Reference  2.  In  the 
original  CM3DT,  the  orientation  of  the  nosetlp  axis  to  the  afterbody  axis  Is 
defined  by  the  angles  <5-|,  <$2,  and  are  inputs  to  the  code,  where  6^ 

and  <S2  define  the  pitch  and  yaw  orientations  of  the  nosetlp  axis  relative  to 
a  bent  afterbody  axis,  and  Is  the  bend  angle  of  the  afterbody  axis.  When 
using  QUICK  with  CM3DT,  due  to  the  plane  of  symmetry  assumption  and  moving 
axis  assumptions  of  QUICK,  <$2  and  6^  must  be  zero  (set  internally  by  CM3DT), 
and  <5^  is  computed  automatically  from  Equation  (2.6). 
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In  using  CM3DT  with  the  QUICK  geometry  definition,  the  angle  of 
attack  (a)  and  sideslip  angle  (0)  are  defined  In  the  Input  relative  to  the 
xq  axis.  With  i1  0  and  <$2  ■  *  0,  the  pitch  and  yaw  angles  relative  to 

the  noset Ip  x  (denoted  by  a^,  0^)  become 

0N  -  0  (2.7) 

*  tan“l  [(sin  ccsa  +  cosfi^  slna)/ 

fcos  6]  cosa  -  Sln61  slna)]  (2.8) 

(Although  both  the  QUICK  and  STEIN  codes  assume  a  pitch  plane  of 
geometric  symmetry  and  0  ■  0°,  the  version  of  CM3DT  modified  to  use  the  QUICK 
geometry  system  requires  only  the  assumption  of  geometric  symmetry.  This  new 
version  of  CM3DT  retains  the  ability  to  treat  non-zero  sideslip.  However, 
when  using  rM3DT  to  generate  Initial  data  for  the  STEIN  code,  geometric  sym¬ 
metry  about  the  pitch  plane  and  0*0°  must  be  assumed.) 

Choosing  x  *  0,  y  «  0,  z  *  0  In  the  nosetlp  coordinates  at  the 

v*  C  C 

point  corresponding  to  xq  «  $b0*  yq  *  0,  zq  »  d-j  results  In  the  following  re¬ 
lationships  between  the  QUICK  and  CM3UT  coordinate  systems: 


xc  a 

(xq 

-  *b0)  cos 

61  -  ^ZQ 

-  d-j )  sin  6^ 

(2.9) 

*c  * 

(Xq 

-  *b0)  sin 

«1  +  (ZQ 

—  d^ )  COS 

(2.10) 

zc  * 

-yQ 

(2.11) 

XQ  * 

*b0 

+  x„  cos 
c 

61  + 

sin  61 

(2.12) 

yq  ■ 

-zc 

(2.13) 

ZQ  * 

dl  - 

x„  sin  0. 

C  1 

+  yc  cos 

51 

(2.14) 

In  the  CM3DT  code  it  is  necessary  to  have  the  body  geometry  defined 
in  the  form  y  *  yfa  (x,  <j>).  To  develop  this  relationship  for  a  geometry  defined 
In  the  QUICK  system  an  Iteration  Is  required,  as  outlined  below. 
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Assuming  that  values  of  x  and  d  are  specified,  an  initial  guess  is 
made  for  y,  and  the  corresponding  values  of  (xcl,  yc,  zc)  are  computed  using 
Equations  (2.3)  -  (2.5).  Values  for  (xq  ,  yq  ,  zq  )  are  then  determined  from 
Equations  (2.12)  -  (2.14).  The  corresponding  values  of  r  and  0  are  in  turn 
determined  from  Equations  (2.1)  -  (2.2). 

It  is  important  to  note  that  CM3DT  and  QUICK  use  opposite  half¬ 
spaces  in  treating  problems  with  a  pitch  plane  of  symmetry.  Thus,  for  values 
of  <p  in  the  range  0  <  <j>  <_  tt  the  values  of  <p  computed  as  outlined  above  will 
fall  outside  of  the  admissible  range  of  0  ,  -tt/2  £0  <_  tt/2.  Since  QUICK 
assumes  a  pitch  plane  of  symmetry,  however,  the  angle  0*  can  be  defined  such 
that  r  (xq,  0)  =  r  (xq,  0*),  where  -tt/2  <_  0*  <_  tt/2.  As  can  be  seen  in  Figure 
2.4,  0*  may  be  defined  as 

-i  ZQ  -  d 

0*  =  tan  -  (2.15) 

Q 

ZQ 
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Once  0*  has  been  determined,  the  actual  body  radius  is  found  from 
the  QUICK  geometry  routines  as  r  =  b  (xq  ,0*). 

This  value  of  b  is  compared  to  the  value  of  r  computed  from  the 
originally  assumed  value  of  y;  if  b  and  r  are  equal,  the  Iteration  is  stopped. 
If  b  and  r  differ,  another  guess  is  made  for  y,  and  the  iterative  process  is 
continued. 


As  described  in  Reference  1,  the  CM3DT  code  makes  use  of  a  trans¬ 
formed  computational  coordinate  system  that  is  fitted  to  the  nosetip  geometry. 

An  iterative  process  is  used  in  CM3DT  to  determine  the  image  of  the  body  sur¬ 
face,  based  on  a  definition  of  the  body  geometry  of  the  form  y  =  yb  (x,  <f>). 

Thus,  when  using  the  QUICK  geometry  option  in  CM3DT,  a  double  iteration  is 
required  to  relate  the  body  surface  in  QUICK  coordinates  to  the  (x,  y,  <f>) 
cylindrical  system  ar.d  then  to  the  transformed  computational  coordinate  system. 
Although  this  double  iteration  is  not  strictly  required  when  using  the  QUICK 
geometry  package  in  CM3DT,  its  use  allows  the  retention  of  the  other  geometric 
definition  options  available  in  CM3DT,  which  are  described  in  Reference  2. 

In  the  definition  of  the  body  geometry  in  CM3DT,  it  is  also  necessary 
to  define  the  body  slopes  9yb/9x  and  3y^/34>.  These  values  can  be  determined  in 
terms  of  the  slopes  9b/9xq  and  9b/90  ,  which  are  provided  from  the  QUICK  ge¬ 
ometry  routines.  By  relating  the  expressions  for  the  outward  body  normal  in 
both  the  CM3DT  and  QUICK  coordinate  systems,  the  desired  body  slopes  are  de¬ 
termined  to  be 

3yb  sin  6^  -  cos  6-j  N1 

__=  |j  (2.16) 

fyb  =  ^  [sin  6-j  sin  <j>  N1  +  cos  4>  N2  + 

90 

cos  5-j  sin0  Nj] 


(2.17) 
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COUPLING  OF  CM3DT  AND  STEIN  FLOW 
FIELD  CODES 

The  STEIN  afterbody  code,  described  in  References  3  and  5,  solves 
the  steady  Inviscid  flow  equations  using  a  forward-marching  (in  space)  explicit 
finite  difference  procedure.  For  this  approach  to  be  applicable,  the  velocity 
component  in  the  marching  direction  must  be  supersonic  at  all  points  in  the 
shock  layer.  In  addition,  this  forward-marching  procedure  requires  the  speci¬ 
fication  of  an  initial  data  surface  which  defines  the  complete  flow  field  at 
a  given  station.  This  section  describes  the  procedure  used  to  y*  .erate  initial 
data  surfaces  for  use  STEIN  afterbody  code  from  CM3DT  nosetip  flow  field  solu¬ 
tions. 


The  STEIN  code  is  formulated  in  a  generalized  coordinate  system 
which  is  fitted  to  body  cross-sections  through  the  use  of  conformal  transfor¬ 
mations.  These  conformal  transformations  are  applied  in  the  STEIN  code  for 
axial  stations  (z  =  constant)  donwstream  of  a  user-specified  point,  z  =  ZMAP2 ; 
no  conformal  transformations  are  used  for  stations  upstream  of  another  user- 
specified  point,  z  =  ZMAP1.  For  ZMAPl<z<  ZMAP2,  the  conformal  transformations 
are  introduced  gradually  to  provide  a  smooth  transition  to  the  complete  mappings. 
When  the  conformal  transformations  are  not  applied  (i.e.,z<  ZMAP1 )  a  polar 
(r,  G)  mesh  is  used  within  each  cross-section,  with  mesh  points  equally  spaced 
in  O  and  equally  spaced  in  r  between  the  body  and  shock  surfaces.* 

For  simplicity,  it  is  assumed  in  this  coupling  effort  that  the  initial 
data  plane  to  be  generated  for  the  STEIN  code  will  be  at  a  station  where  the 
conformal  mapping  procedure  is  not  used,  i.e.,  at z <  ZMAP1.  Thus,  the  initial 
data  to  be  generated  from  the  CM3DT  solutions  is  required  to  be  in  a  polar 
coordinate  system  within  the  specified  starting  plane. 


*  The  axial  coordinate  z  in  the  STEIN  coordinate  system  is  equivalent 
to  xq  in  the  QUICK  coordinate  system, 
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Inputs  required  to  the  routine  in  which  the  initial  data  surface  is 
generated  are  ZSTART,  NCI,  and  MCI,  where  ZSTART  is  the  initial  value  of  2  for 
the  STEIN  calculation,  NCI  is  the  number  of  points  in  the  radial  direction  and 
MCI  is  the  number  of  points  in  the  circumferential  direction.  (These  parameters 
have  the  same  definition  in  this  application  as  in  their  use  as  input  para¬ 
meters  to  the  STEIN  code.)  With  these  parameters,  the  STEIN  mesh  points  are 
located  in  polar  coordinates  at  z  *  ZSTART  as 

©m  a  -  tt/2  +  (m-l)/(MCl-  1 )  it  *,  (3.1) 

m  =  1 ,  2.  .  .  .,  MCI 

rn,  m  *  b(z,  Om;  +  (n  -  1 )/ (NCI  -  1)  [c(z,0m)-b(z,0m)J 

n  =  1,2,  ....  NCI  (3.2) 


where  r  =  b  (2,  0)  represents  the  body  surface  and  r  =  c  (z,  0)  represents 
the  bow  shock  surface. 

The  first  step  in  the  generation  of  the  Initial  data  is  to  determine 

the  location  and  slopes  of  the  body  points  In  the  initial  data  plane  from  the 

QUICK  geometry  routine.  Next,  it  is  necessary  to  locate  the  bow  shock  points 

r  =  c  (z,  0);  since  the  bow  shock  location  is  known  in  CM3DT  coordinates  as 
m 

y  =  ys  (x,  <)>),  an  iterative  interpolation  procedure  is  required.  The  rela¬ 
tionships  between  the  STEIN  (QUICK)  and  CM3DT  coordinate  systems  used  In  this 
procedure  are  detailed  in  Section  2. 


Once  the  bow  shock  point  r  =  c  (z,  0m)  has  been  located  in  CM3DT 
coordinates,  the  bow  shock  slopes  3y$/3x  and  3ys/3<{>  may  be  readily  deter¬ 
mined  through  interpolation.  These  shock  slopes  may  be  expressed  In  the  STEIN 
coordinate  system  as  3c/9z  and  3c/30  ,  which  can  be  computed  from 

3c/3z  =  -  N2/N1  -  dz  (sin  0  -  3c/30  )  (3.3) 

3c/30  =  -  c  N3/N1  (3.4) 

where 


-  cos©  sin  <p  +  cosG  cos<(i  ys  /y$  + 
sin©  sin  6-|  ySx  +  sin©  cos  6^  (cos<J>  + 
sin  $  yS(p  /ys) 


(3.5) 
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(3.6) 


N2  *  -  y$x  cos  <5^  +  sin  5^  (cos  4>  +  sin  4>  yS({/ys) 

*  sin©  sin  4>  -  sin©  cos$  ys^/y$  + 
cos©  sin  6j  y$x  +  cos©  cos  (3.7) 

(cos  <t>  +  sin  4  ys  /y$) 

In  using  the  above  formulae,  care  must  be  taken  to  ensure  that  the 
correct  values  of  4>  and  3ys/34>  are  used.  For  a  given  point  (r,  0,  z),  where 
-  tt/2  £  ©  <_  tt/2,  the  corresponding  point  In  the  (x,  y,  4>)  CM3DT  coordinates 
will  lie  in  the  range  it  <_  4)  <_  2n.  The  Interpolation  process  Is  carried  out 
in  the  CM3DT  system  using  the  circumferential  angle  <J>*»  where,  with  a  pitch 
plane  of  symmetry, 


<J>*  a  2  TT  -  <p 

(3.8) 

<_ tt.  From  the  symmetry  conditions. 

ys  (x,  4)  «  y$  (x,  <<>*) 

(3.9) 

ySx  (x.  4»)  s  ySx  (x,  ♦*) 

(3.10) 

ys^  (x.  ♦)  *  -  ys<p  (x,  ♦*) 

(3.11) 

Once  the  body  and  bow  shock  points  have  been  located  in  the  STEIN 
initial  data  plane,  the  location  of  all  STEIN  mesh  points  in  this  plane  can 
be  easily  determined  from  Equations  (3.1)  «ind  (3.2).  Knowing  (r,  0,  z)  in 
the  STEIN  coordinate  system,  the  corresponding  (x,  y,  4>)  in  CM3DT  coordinates 
is  easily  determined  using  the  transformations  presented  In  Section  2.  With 
the  corresponding  (x,  y,  4>)  known  for  a  given  mesh  point  In  the  Initial  data 
plane,  the  dependent  variables  of  the  CM30T  flow  field  solution  at  that  point 
are  readily  determined  by  interpolation. 

As  described  in  Reference  1,  the  CM3DT  flow  field  variables  are 
P,  u,  v,  w,  and  s/R,  where  P  *  log  p  and  u,  v,  w  are  the  velocity  components 
in  the  (£,  n,  ©)  transformed  space.  Prior  to  interpolating  for  the  required 
flow  field  data,  the  CM3DT  velocity  components  are  transformed  to  the  (x,  y,  4>) 
velocity  components,  written  as  U,  V,  W,  respectively,  as 
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U  «  uC  -  vS  (3.12) 

V  *  uS  +  vC  (3.13) 

W  «  w  (3.14) 

where  C  +  IS  s  e"1u,  as  defined  In  Reference  1. 

It  must  again  be  noted  that  the  Interpolation  Is  carried  out  In 
terms  of  0*  *  2 tt  -  0,  where  0  £  0*  <_  tt  .  The  symmetry  conditions  applied  to 
the  CM3DT  flow  variables  are 

P  (x,  y,  0)  ■  P  (x,  y,  0*)  (3.15) 

s/R  (x,  y,  4>)  *  s/R  (x,  y,  0*)  (3.16) 

U  (x,  y.  ♦)  -  U  (x,  y.  0*)  (3.17) 

V  (x,  y,  <J>)  -  V  (x,  y,  0*)  (3.18) 

W  (x,  y,  0)  *  -W  (x,  y,  0*)  (3.19) 

In  STEIN  the  dependent  thermodynamic  variables  used  in  the  calculation 
are  defined  as 

P  =  log  (p/p^  )  (3.20) 

s  *  (5  -  sJ/Cv.  (3.21) 

where  (  )  denotes  a  dimensional  quantity.  The  STEIN  velocity  components  are 

defined  as  us>  vs,  and  ws,  which  correspond  to  the  Cartesian  directions  yq  , 

Zq,  xq  (in  QUICK  coordinates),  respecti vely.  The  dimensional  velocity  com¬ 
ponents  are  defined  from 

u$  *  -  V  sin0  -  W  cos  0  (3.22) 

vs  *  -  U  sir6-|  +  V  cos6^  cos0  -  W  cos6-|  sin0  (3.23) 

ws  8  U  cos5^  +  V  s1n<5^  cos0  -  W  s1n6-j  sin0  (3.24) 

and  the  non-dimensional  values  are  determined  from 

us*  V  ws  *  V  V  V'/p7£  <3-25) 


16. 


In  determining  the  starting  line  data  the  quantities  P,  s,  us,  and 
vs  are  determined  by  the  Interpolation  procedure.  The  remaining  velocity  com¬ 
ponent,  ws,  is  determined  from  the  requirement  of  constant  total  enthalpy  in 
a  steady  inviscid  flow,  from 

ws*[  2  (Hq  -  h)  -  (us2  +  vs2B1/2  (3.26) 

where  HQ  is  the  total  enthalpy  and  the  static  enthalpy  is  obtained  in  the  form 
t  ■  h  (p,  s). 

After  determining  the  required  initial  data  surface  for  the  STEIN 
code  from  the  CM3DT  solution.  It  Is  necessary  to  evaluate  the  force  and  moment 
integrals  on  the  nosetip,  and  relate  these  to  the  STEIN  coordinate  system  if  the 
subsequent  STEIN  calculation  is  to  Include  force  and  moment  computations.  The 
integrations  are  carried  out  In  the  CM3DT  coordinate  system,  and  It  is  thus 
necessary  to  define  the  downstream  limits  of  Integration  of  the  form  x  «  (d>) , 

which  is  determined  by  transforming  the  body  points  in  the  STEIN  initial  plane 
to  (x,  y,  4>)  coordinates  and  interpolating  to  find  the  values  x£  in  the  <p  planes 
used  in  the  CM3DT  calculation. 

With  the  assumptior  of  a  pitch  plane  of  symmetry  only  two  forces  are 
present,  the  axial  and  norma'  forces.  These  forces  may  be  evaluated  as  (assuming 
zero  base  pressure) 

?AN  -  2  >1  /JeM  p  yb  y„x  d*  dp  (3.27) 

F"N  '  '2  ’l  i&'p  yfa  (cost  +  sint  /yb  )  dx  dt  (3.28) 

The  pitching  moment  (which  is  the  only  non-zero  moment  with  a  pitch  plane  of 
symmetry)  is  referenced  to  a  point  (xcp  yep  zq  )  in  the  CM3DT  Cartesian  co¬ 
ordinate  system.  The  point  (xCp  yep  zq )  is  determined  from  the  transforma¬ 
tion  from  QUICK  coordinates  of  the  moment  reference  point  (xQp  yQp  zq-j),  which 
are  defined  in  the  STEIN  code  through  the  inputs  VMO  (3)  =  xq ,  VMO  (1)  «  yQ^ , 
and  VMO  (2)  3  zq^ .  In  STEIN  it  is  assumed  that  VMO  (1)  =  yq^  3  0;  thus  it 
follows  that  Zq  *  0.  The  resulting  expression  for  nosetip  pitching  moment  is 
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M  ■  2  /q  pyb  [  X  (cos<f>  +  Sin4»/yb  yb^  ) 


+  yb  ybx  cos  <j>]  dx  dd>  (3.29) 

+  Xc1  FN  “  ycl  fa 

Finally,  the  axial  and  normal  forces  in  CM3DT  coordinates  can  be 
related  to  the  STEIN  coordinate  system  through 

Fa  *  co$61  Fan  +  s1n61  Fnn  (3.30) 

Fjj  *  -  sin6j  F/^  +  cosS^  Fnn  (3.31) 

In  addition  to  the  nosetlp  force  and  moment  integrals,  the  initial 
data  for  STEIN  also  includes  the  nosetip  surface  area.  This  quantity  Is  deter¬ 
mined  from  evaluating  the  Integral 

A  *  2  f  /XE(*)  yb  /TT^  "+'yb//yb2"  dx  d4>  (3.32) 

The  initial  data  generated  for  the  STEIN  code  from  a  CM3DT  solution 
is  written  on  a  binary  file  in  a  for  't  compatible  with  the  existing  starting 
procedure  in  the  STEIN  code.  The  processes  required  to  start  a  STEIN  calcula¬ 
tion  using  this  procedure  are  described  in  the  following  section. 
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SECTION  4 


USER'S  GUIDE  TO  THE  MODIFIED  CM3DT  CODE 


The  process  of  making  a  STEIN  calculation  with  CM3DT-generated 
initial  data  requires  three  steps.  First,  the  geometry  of  Interest  Is  de¬ 
fined  using  program  QUiCK,  which  produces  an  Intermediate  file  or  deck  of 
geometry  data.  Second,  the  modified  version  of  CM3DT  (which  includes  the 
subset  of  QUICK  routines  known  as  SUB-QUICK)  reads  the  Intermediate  geometry 
file,  performs  the  nosetlp  calculation,  and  generates  an  output  file  of 
initial  data  for  STEIN.  Finally  the  STEIN  program  proceeds  from  the  initial 
data  plane  generated  by  CM3DT.  QUICK  and  STEIN  are  documented  in  References 
4  and  5 ; thi s  section  describes  modifications  to  CM3DT  input  and  operation. 


Load  Sequence 


CM3DT,  when  combined  with  SUB-QUICK,  occupies  250000  octal  words 
of  memory  on  a  CDC  Cyber  176.  This  load  sequence  includes  subroutine  GC0M6 
which  generates  geometry  using  the  QUICK  intermediate  data  file,  but  not 
subroutines  GE0M3,  GE0M4,  GE0M5,  or  ACOEF,  which  treat  other  geometry  options 
in  CM3HT  and  art  not  used  in  this  application.  Also  added  to  the  load  se¬ 
quence  art  iuoroutlnes  STSTEIN  and  CLCINT,  which  generate  STEIN  starting 
line  deta  and  perform  a  general  purpose  integration,  respectively.  In  addi¬ 
tion,  toe  estimate  of  memory  allocation  presumes  the  loaoing  of  either  sub¬ 
routine  PVE3DL  or  subroutine  PVE3DNC  but  not  both,  as  discussed  in  Reference 
2.  Thui,  the  list,  cf  subroutines  to  be  compiled  for  execution  of  the  modi¬ 


fied  CM3DT  program  is  as  follows: 


CM3DT 

PVE3DL  or  PVE3DNC 

HING3D 

GRID3D 

0UT3D 

BFL3CM 

MAP  3D 

GE0M6 

STSTEIN 

CLCINT 


SHOCK  \ 
STAUPT  ' 
RGAS 
CUFT1 
TBL1  ' 
IDEAL 
RLERR 
TLU1 
SHKTAB 
FREE 
ATMP 
CUFT2 
ATERR 
LINB 
TBL  , 
AIR  / 
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>  CM3DT 


> 


GEOMIN 
CSGEOM 
CURVES 
BLMSET 
CSMSET 
BLGEOM 
VDOTV 
MDOTV 
THELIM 
CSMCOE 
CSMINT  ! 
CSCALC 
CSMFLT  j 


SUB-QUICK 


Modifications  to  CM3DT  Input 


Modifications  to  the  CM3DT  Input  consist  of  two  types.  Those 
changes  which  are  necessary  for  the  Inclusion  of  the  QUICK  geometry  capa¬ 
bility  are  incorporated  Into  namelist  INPUT,  while  additional  variables 
necessary  for  the  generation  of  initial  data  for  STEIN  comprise  a  new  name- 
list  ST,  which  follows  INPUT.  Definitions  of  input  variables  which  are  new 
or  have  been  altered  are  listed  below: 


Namelist  $  INPUT 

ALPHA  Angle  of  attack  (degrees)  relative 

to  QUICK  axis 

BETA  Sideslip  angle  (degrees)  relative 

to  QUICK  axis 

IQUICK  f  0  when  QUICK  geometry  used 

(default  0) 

XBB(L)  Xg  values  at  4>  *  PHIBD(L )  defining 
downstream  boundary  (NPHI  values; 
PHIBD's  generated  automatically 
for  NPHI<  3  or  NPHI  *  4  and 


XBO 

XHO 


ISA  «  2) 

Xg  value  of  most  forward  nosetip  point 
Xg  value  used  to  define  point  on  com¬ 
putational  axis  downstream  of  XBO 


N.B.:  When  using  CM3DT  to  generate  initial  data  for  STEIN,  must  have 
ISA  *  1  (default)  with  BETA  -  DEL 2  »  DELA  -  0.0.  DELI  will  be 
defined  internally  by  CM3DT. 
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Namelist  $ST 


ZSTART  Axial  location  (In  QUICK  coordinates)  of  Initial 
data  plane 

MCI  Number  of  points  to  be  generated  In  the  circumferen¬ 

tial  direction 

NCI  Number  of  points  to  be  generated  In  the  radial 

direction 

IAERO  «  0  no  aerodynamic  coefficients  to  be  computed  In 
STEIN 

t  0  generate  Initial  force  and  moment  values  for 
STEIN 

XQREF  Axial  location  of  moment  center  In  QUICK  coordi¬ 

nates  (equal  to  VM0(3)  In  STEIN  Input) 

ZQREF  Vertical  location  of  moment  center  In  QUICK  coordi¬ 
nates  (equal  to  VM0(2)  In  STEIN  Input) 

XQ1.D1XQ1  Axial  and  radial  locations  of  the  geometric  stagna¬ 
tion  point  In  QUICK  coordinates.  (These  values 
need  be  input  only  for  a  restarted  CM3DT  solution, 
and  can  be  obtained  from  the  section  of  CM3DT  out¬ 
put  displaying  QUICK  geometry  parameters.) 


File  Structure 


In  addition  to  the  CM3DT  solution  file  (TAPE21)  two  additional 
files  must  be  declared  local  to  the  modified  CM3DT  program:  the  QUICK 
gemoetry  data  file  (TAPE22)  and  the  STEIN  initial  data  file  (TAPE23).  The 
existence  of  the  CM3DT  solution  file  enables  the  generation  of  STEIN  ini¬ 
tial  data  without  repeating  the  entire  calculation.  This  is  done  by  re¬ 
starting  the  CM3DT  solution  and  setting  KSTART  equal  to  KMAX.  By  this 
action,  control  is  passed  immediately  to  the  initial  data-generating  sub¬ 
routine  STSTEIN.  It  should  be  noted  that  implementation  of  this  pro¬ 
cedure  requires  that  the  user  include  as  input  the  variables  XQ1  and  D1XQ1, 
since  these  are  not  written  on  the  solution  file. 
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Sample  Cases 


Two  sample  cases  are  listed  below.  The  first  Is  a  CM3DT  calcu¬ 
lation  at  Mach  24  and  20°  angle  of  attack  of  the  Shuttle  Orblter  with 
generation  of  starting  line  data  for  STEIN  at  *  •  44.44.  The  second  case 
is  similar  to  the  first,  except  that  the  STEIN  initial  data  Is  generated 
from  an  already  completed  CM3DT  solution. 


Sample  Case  1 


JOB  • _ _ 

ACCOUNT* _ _ 

REQUEST *  TAPES  1  * ♦PF. 

REQUEST*  TRPE£3* ♦PF. 

ATTACH* OLD 1*CM3FL* ID* _ 

UPDATE* P-OLD1. 

FTN<I*B=CM3DT*L*0> 

ATTACH* OLD£* QU I CKPL* ID* _ 

UPDATE* P=OLD£. 

FTN<I*B*QUICK*L«0> 

AT T ACH *  TAPE££ *  QU I CRD AT A *  ID* _ 

LDSET (PRESET *ZERO> 

L0RD<CM3DT*QUICK> 

EXECUTE. 

EXIT*  U. 

CATALOG* TAPE£1* CMS* ID* _ 

CATALOG*  TAPE£3*  CN3ST* ID- _ 

/  §  /  Q 

♦COMPILE  CM3DT *  PVE3DL *  H I NG3D* GR I  DSD *  MAP3D *  GE0M6 
♦COMPILE  BFL3CM* 0UT3D*  STSTEIN*  CLCINT*  STAGPT*  SHOCK 
♦COMPILE  RGAS*  CUFT1 * TBL1 *  IDEAL* RLERR*  TLU1 *  SHKTAB 
♦COMPILE  FREE* RTMP* CUFT£* ATERR* LIMB* TBL 

♦COMP I LE  CURVES » GEOM IN*  CS6E0M*  BLMSET * CSMSET *  BL6E0M *  VDOTV 
♦COMP I LE  MDOTV *  THEL I M  *  CSMCOE *  CSM I NT  *  CSCALC » CSMFLT 

7  /0/C| 

CM3DT  <LAMBDA>  SHUTTLE  W-  QUICK  GEOMETRY  MACH  £4  ALPHA*£0 

$ INPUT 

NMRX*E  *  MM AX* 1 8  *  LMAX*? * 

KMAX*500* KOUT=500* KTAPE=£50* 

AMI NF*£4.  *ALPHA=£(i.  * 

NPH I *£  *  XBB*6 0 . *60. * XB 0= 0 . *  XH 0*15 0 . * 

RN=50. * JA=7* I LAM* 1  * IQUICK*1$ 

tST 

ZSTART=44. 44*  MCI *31  *  NCI *15* 

I  AERO*  1  *  XG!REF=840.  7*  ZQREF=375.  $ 
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Sample  Case  2 


JOB* _ _ 

ACCOUNT, _ 

ATTACH. TAPES9, CM3, ID* _ 

COPY  <TAPE£9, TAPES 1) 

RETURN <TAPE£9) 

REMIND  (TAPE£1) 

REQUEST , TAPES3, ♦PF. 

ATTACH, OLD 1 , CM3PL, ID* _ 

UPDATE, P-OLDl. 

FTN  (I , B*CM3DT , L=0) 

ATTACH, OLDS, QUICKPL, ID* _ 

UPDATE, P=OLD£. 

FTN  <1 , B*QUICK, L=0) 

ATTACH, TAPE££, QUICKDATA, ID* _ 

LDSET  (PRESET=ZERO) 

LOAD (CM3DT, QUICK) 

EXECUTE. 

EXIT, U. 

CATALOG,  TAPES3,  CM3ST ,  ID* _ _ 

?'8/9 

♦COMP I LE  CM3DT , PVE3DL >  HING3D, 6RI DSD , MAP3D , GEOM6 
♦COMPILE  BFL3CM, 0UT3D, STSTEIN,CLCINT, STAGPT, SHOCK 
♦COMPILE  RGAS, CUFT1 , TBL1 , IDEAL, RLERR, TLU1 , SHKTAB 
♦COMPILE  FREE, ATMP, CUFTS, ATERR,LINB, TBL 

7/8/9 

♦COMPILE  CURVES, GEOMIN, CSGEOM, BLMSET , CSMSET , BLGEOM, VDOTV 
♦COMP  I LE  MDQTV, THEL I M, CSMCOE , CSM I NT , CSCALC , CSMFLT 
7/8/9 

CM3DT  (LAMBDA)  SHUTTLE  M/QUICK  GEOMETRY  MACH  £4  ALPHA*£Ci 

$  INPUT 

NMAX=6, MMAX=18, LMAX=9, 

KM AX =5 0  0 , KSTART =5 0 0, 

AMI NF=c’4 .  ,ALPHA=£0.  , 

NPHI*£,XBB=60. ,60.  ,XB0*0.  ,XH0=150. , 

RN*50. , JA*7 , ILAM*1 , IGUICK=li 
f.ST 

2START=44 . 44 , MC 1 *3 1 , NC 1  *  1 5 , 

XQ 1=0. , D 1 XQ 1 =335 . , 

I  AERO* 1 , XQREF-840. ?,ZQREF=3?5. * 

6/ 7/8/9 
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